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RiiLATIONSHIPS  B^TiVEEN  THE  WIND-FIELD,  THE 
TRAl'JSPORT  FIELD  AND  THE  MEAN  MASS-FIELD 
IN  A  NON-HOMOGENEOUS  OCEAN 


V.B. SHTOKMAN 


(Presented  by  P.P.SHIRSHOV, 

Member  of  the  Academy,  December  9,  1947) 


The  classical  theory  of  the  horizontal  circulation  in  a 
hcanogeneous  ocean,  as  constructed  by  V.EKMAN  {l},  suffers  from  the 
defect  that  it  considers  only  the  frictional  forces  acting  in  hori¬ 
zontal  planes.  In  consequence,  once  the  water  on  the  sea-bottom 
has  acquired  a  motion,  vertical  velocity-gradients  must  be  developed 
just  as  great  as  those  in  the  surface ‘layer  of  the  sea.  This  result 
does  not  agree  with  observed  data,  which  show  that  the  vertical 
velocity-gradients  in  a  non-homogeneous  ocean  are  rapidly  damped  out 
with  increasing  depth,  and  at  the  bottom  are  to  all  practical  in¬ 
tents  and  purposes  insignificant.  This  rapid  attenuation  of  the 
vertical  velocity-gradients  is  a  result  of  intensive  horizontal  ex¬ 
change  and  transfer  of  motions,  setting  up,  in  vertical  planes 
parallel  to  the  horizontal  components  of  flow,  large  frictional 
forces  which  almost  entirely  counter-balance  the  tangential  stresses 
due  to  the  wind  on  the  surface  of  the  ocean.  Thus  the  distinctive 
feature  of  an  established  motion  in  a  non-homogeneous  ocean  is  the 
fact,  confirmed  by  observation,  that  the  frictional  forces  at  the 
bottom  are  extremely  small  in  comparison  with  the  stresses  set  up 
by  the  turbulent  "side-friction”.  On  the  basis  of  these  findings, 
the  equations  for  the  components  of  the  steady— state  transport  • 
in  a  non-homogeneous  ocean  may  be  written  in  the  following  fom: 


A: 


-j-Tj,  —  cpSj,=g 


dy  ’ 


(1) 


*  "Transport"  in  this 
(quantity  of  flow). 


translation  corresponds 
as  defined  in  equations 


to  Ekman’s  Strommenge 
la,  (Translator) 


Here  Ai  represents  the  coefficient  of  horizontal  turbulent  friction 
(which  is  of  the  order  of  lOP  CGS),  o  =  Soo  sin  <p  is  the  Coriolis 
parameter,  g  is  the  acceleration  of  gravity,  Tx  and  Ty  are  the  hori¬ 
zontal  and  y-components  of  the  tsingential  wind-pressure  at  the 
surface  of  the  ocean,  p  is  the  mean  density  of  the  water  in  a  column 
extending  from  the  surface  (2=0)  to  the  bottom  (s  -  h),  and  Sx  and 
3y  are  the  components  of  the  transport; 

h  h 

Sj,—  ^uc/z,  Sy  —  ^vdz,  (la) 

0  0 


where  u  and  v  are  the  horizontal  coniponents  of  the  velocity  of  flow. 
Again,  we  have 


The  equation  of  continuity  for  the  transport 
may  be  written  in  the  form: 


dx 


(.pS,)+-^(pS,)  =  0. 


Since  the  variations  of  p  along  x  and  y  are  very  small  in  comparison 
with  the  corresponding  variations  of  Sx  and  3y,  we  may  neglect  the 
quantity  p  in  the  last  equation,  and  write 


For  the  same  reason  we  may,  in  differentiating  equations 
(1)  with  respect  to  x  and  y,  put  p  approximately  equal  to  xmity« 
Eliminating  the  quantities  dQ/dx  and  dQ/dy  from  equations  (l),  we 
have  by  (s): - 


v2rot5=-rot7’//l;,  (3j 

where 

V^s=d^/dx^+d-/dy\  roXT^dTJdx  —  dTJdy,  io\S=dSJdx-dSJdy. 

If  we  introduce  the  transport  function  ijf,  related  to 
and  components  of  the  transport  by  the  expressions:- — 
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(4) 


5^  =  —  d^/dy,  Sy  =  d^/dx, 

then  (3)  is  transformed  into  the  equation 


y-»v|;=—  rOtT/yAj, 


(5) 


analogous  to  the  equilihriura-equations  of  a  plate  under  the  action 
of  applied  loads  (to  which  the  quantity  rot  T  is  analogous. )  In 
equation  (5),  v*  represents  the  biharmonic  operator.  The  boundary 
conditions  for  the  transport  on  the  contour  r  of  the  ocean  shore  are 
analogous  to  those  for  a  plate  held  by  pinching  it  edge-on: - 


where  n  is  the  normal  to  the  contour.  Conditions  (6)  are  derived  by 
considering  that  the  tangential  and  normal  components  of  the  trans¬ 
port  at  the  shore-line  must  reduce  to  zero.  Let  us  now  find  the 
equation  relating  the  mass-field  to  the  wind-field.  Finding  S^.  from 
the  first  equation  of  (l)  and  substituting  it  in  the  second  equation, 
we  solve  the  latter  for  Syi 


dx 


dy 


c 


Similarly: 


cS, 


g  47  +ry+  iJC  ^  v'Q  -  ^  vV,-  v*5^. 

oy  c  ax  c  c 


(B) 


Di^’ferentiating  equation  (7)  with  respect  to  y,  equation 
(8)  v/ith  respect  to  x,  and  combining  the  results,  we  get: 


fCliv5=rot7'-f  y^Q 


c 


div  7* - y*div5 

c 


or,  by  condition  (2): - 


J 


(9) 


Thus  the  dis trihut ion  of  masses,  represented  hy  Q,  de¬ 
pends  not  only  on  the  Coriolis  parameter  and  rot  T,  hut  also  on  the 
quantity  v^div  T  =  div  grad  div  T,  The  latter  is  usually  very  small 
in  comparison  with  the  first  term  of  (9),  However,  it  cannot  he 
neglected  in  regions  of  marked  convergence  or  divergence  of  air- 
currents  (for  instance  the  center  of  a  cyclone  or  anti-cyclone). 


If  to  describe  the  mass-field  we  ^hoose  to7  use 
the  quantity  P,  related  to  the  dynamic  height  hy  the 
expression 

h 

P=l  Ddz 

0 

where  the  ^^eopotentia]J7  height  JD  is  measured  in  dy¬ 
namic  centimeters,  then  the  last  equation  will  he 
written  in  the  form; 


y»Z)  =  _I_  ,  div  7’')  .  (1 0) 

lO'A  ly.  ^  j 


The  boundary  conditions  to  which  Q  or  D  are  subject  will 
he  found  from  equations  (7)  and  (8),  If  at  a  given  point  we  take 
X  as  the  direction  tangential  to  the  shore-line  and  y  as  the  normal 
thereto,  and  if  we  hear  in  mind  that  at  the  shore-line  =  Sy  =  0, 
then  we  obtain  the  following  conditions  for  Q: 


dll  T72  ^ 


dx  ~  '  c 


dy 

dQ 


dy  '  c  dx  c 
Recalling  (4),  we  may  write; 


r 

c  ^ 


(11) 

(12) 


* 


-  4  - 


Replacing  by  its  expression  in  terras  of  rot  T  as 
given  by  (5),  we  obtain 


_d_ 

dx 


rot  7; 


aI  a 


dy 


rot  7. 


(13) 


Substituting  (13)  in  (11)  and  (12),  we  obtain  equivalent 
boundary  conditions  in  the  following  form; - 


_ Zk _ -dl  ~ 

dx  g  c  dy  “  eg  ^  c  dx 


(14) 


;ii— rotr. 

dy  g  '  e  dx  eg  eg  dy 


At  .  ao 


(15) 


the  V2  sign  in  (14): 


Finally,  let  us  substitute  expression  (15)  for  under 


dx  g 


— 

c 


+  ^  V 
\  g  c 

+  4rv‘‘7V- 


iiQ 

dx 


At 


eg 


At  ^Z'r  _  AkA.  rot7. 


eg 


eg  dx 


Expanding  the  operators  and  reducing  similar  terms,  we  see 
that  all  terms  except  the  first  in  the  right-hand  side  of  this 
equation  cancel  each  other#  Thus  we  obtain  a  very  simple  boundary 
condition  on  the  shore-line  contour: 


dQ/dx=rjg  (16) 

and  in  the  same  way,  a  second  condition: - 

dQldy^T,lg.  (17) 


If  we  use  Df  the  boundary  conditions  will  obviously 
be  written  in  the  form; 


A2-=Ll  ■  (18) 

dx  ]0*h  ’  dy  lO’/i  ’ 
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Let  us  now  establish  the  relation  between  the  field  of  masses  and 
the  transport.  Differentiating  the  first  equation  of  (l)  with 
respect  to  the  second  with  respect  to  y,  and  combining  the  re¬ 
sults  (p  =  l),  we  obtain 


gSJ^Q=c  rots  +  div  T  +  i4/V*div  S 


or,  by  condition  (8):— 


y2Q__L  (c rots 4-div ri. 
g 


(19) 


Prom  (19)  it  follows  that  the  adaptation  of  the  mass-field 
to  the  transport  field  takes  place  independently  of  the  size  of  the 
coeii^icient  of  horizontal  turbulent  exchange  Ai.  At  the  equator, 
where  o  =  0,  the  adaptation  of  masses  is  regulated  exclusively  by 
the  quantity  div  T,  V/e  note  that  formulae  (7)  and  (8)  may  be  used 
for  calculating  the  transport  components  from  observed  mass-field 
and  wind-field  data,  if  we  substitute  expression  (13)  in  (7)  and  (8). 
In  this  way  we  obtain  the  formulae;  — 


<:  =1.  T'x  I  Mg 

^  c  dx  c  c* 

=  —  s  \ 

'  c  dy  c 


^  a7 


(20) 

(21) 


replacing  the  well-known  formulae  of  Skraan  [2]  and  Jakhelln  {3].  An 
estimate  of  the  order  of  magnitude  of  the  terms  in  (20 )  and  (21)  v/ill 
show  that  the  first  three  terms  on  the  right-hand  side  are  often 
identical  in  amount. 

Thus  if  we  take  "side-friction"  and  wind  into  account,  it 
may  introduce  substantial  corrections  in  Ekman’s  formulae.  When 
there  is  no  v/ind,  or  v/hen  its  effect  may  be  neglected,  expressions 
(20 )  and  (21 )  are  a  good  deal  simplified. 
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